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channels under per-modem power constraints and 
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Vincent Le Nir, Marc Moonen, Jan Verlinden, Mamoun Guenach 

Abstract — In xDSL systems, crosstalk can be separated into two 
categories, namely in-domain crosstalk and out-of-domain crosstalk. In- 
domain crosstalk is also refered to as self crosstalk. Out-of-domain 
crosstalk is crosstalk originating from outside the multi-pair system 
and is also denoted as external noise (alien crosstalk, radio frequency 
interference,...). While self crosstalk in itself can easily be canceled by 
a linear detector like the ZF detector, the presence of external noise 
requires a more advanced processing. Coordination between transmitters 
and receivers enables the self crosstalk and the external noise to be 
mitigated using MIMO signal processing, usually by means of a whitening 
filter and SVD. In this paper, we investigate the problem of finding the 
optimal power allocation in MIMO xDSL systems in the presence of 
self crosstalk and external noise. Optimal Tx/Rx structures and power 
allocation algorithms will be devised under practical limitations from 
xDSL systems, namely per-modem total power constraints and/or spectral 
mask constraints, leading to a generalized SVD-based transmission. 
Simulation results are given for bonded VDSL2 systems with external 
noise coming from ADSL2+ or VDSL2 disturbing lines, along with a 
comparison between algorithms with one-sided signal coordination either 
only at the transmit side or the receive side. 

Index Terms — MIMO systems, Optimization methods 



I. INTRODUCTION 

The growing demand for high speed services in access networks 
calls for new paradigms offering an increased capacity and better 
performance. Thanks to the success of x-Digital Subscriber Lines 
(xDSL) and ADSL in particular, service providers begin to bind 
copper pairs, allowing customers to be served with higher bitrates 
through the usage of adequate Multiple Input Multiple Output 
(MIMO) signal processing algorithms. This processing is also to 
provide a suitable crosstalk interference mitigation. In xDSL systems, 
crosstalk can be separated into two categories, namely in-domain 
crosstalk and out-of-domain crosstalk. In-domain crosstalk is also 
refered to as self crosstalk. Out-of-domain crosstalk is crosstalk 
originating from outside the multi-pair system and is also denoted 
as external noise (alien crosstalk, radio frequency interference,...). 

Self crosstalk cancellation has been studied for two-sided coor- 
dination vector channels and for one-sided coordination Multiple 
Access Channels (MAC) or Broadcast Channels (BC). For two- 
sided coordination vector channels, the Channel State Information 
(CSI) is available at both the transmitter and the receiver. For this 
full vectoring problem, the optimal precoding at the transmitter 
and equalization at the receiver as well as optimal Power Spectral 
Densities (PSD's) are obtained through the SVD of the channel matrix 
combined with standard waterfilling [1]. For one-sided coordination 
MAC or BC, only the receivers for the MAC or the transmitters for 
the BC can cooperate. It has been shown that the optimal structure 
for MAC is a Minimum Mean Square Error-Decision Feedback 
Equalizer (MMSE-DFE) along with a power allocation found by 
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exhaustive search [2]. For BC, a similar (dual) optimal structure has 
also been desribed called MMSE-Dirty Paper Coding (MMSE-DPC). 
Moreover, owing to the diagonal dominance structure of the channel 
matrix, all these optimal structures can be simplified to Zero Forcing 
(ZF) solution for the MAC (or a simple Diagonalizing Precoder (DP) 
for the BC) with transmit PSD's obtained by single-user waterfilling 
[3], [4]. Finally, when there is no coordination, neither at the receive 
side nor at the transmit side, this leads to Interference Channels (IC) 
where spectral management is employed to reduce crosstalk. For IC, 
the optimal transmit PSD's have been found by means of Optimal 
Spectrum Balancing (OSB) [5], [6]. 

External noise is most often the predominant interferer and orig- 
inates from outside the multi-pair system. With external noise, the 
diagonal dominance structure of the channel matrix is destroyed by 
the necessary whitening, and hence the simple ZF solution along 
with single-user waterfilling is found to be suboptimal. The basic 
idea of external noise cancellation is then to exploit the correlation of 
the noise to improve the performance of the transmission and hence 
to increase the total capacity. This noise correlation can appear in 
the spatial domain (between pairs), the frequency domain (between 
tones) or the mode domain (between common-mode and differential- 
mode) [7]. In a recent paper [8], it was shown that there is more 
benefit in exploiting the noise correlation between pairs rather than 
the correlation between tones. 

In this paper, we investigate the problem of finding the optimal 
power allocation in MIMO xDSL systems under self crosstalk and 
external noise and with two-sided coordination exploiting the noise 
correlation between pairs. Coordination between transmitters and 
receivers enables the self crosstalk and the external noise to be 
mitigated using MIMO signal processing by means of a whitening 
filter and Singular Value Decomposition (SVD). Optimal transmit- 
ter/receiver (Tx/Rx) structures and power allocation algorithms will 
be devised with practical limitations from xDSL systems, namely 
per-modem total power constraints and spectral mask constraints, 
leading to a generalized SVD-based transmission. Compared to [8] 
where the external noise is mitigated under a total power constraint 
without self crosstalk, the proposed algorithms mitigate the external 
noise under per-modem total power constraints and spectral mask 
constraints considering the self crosstalk. Capitalizing on the results 
of [8] where it is shown that there is more correlation in the spatial 
domain than the frequency domain, we consider the correlation only 
in the spatial domain. Contrary to [7] where CSI is available only 
at the receive side and the external noise is mitigated by common- 
mode exploitation, our algorithms assume CSI at the transmit and the 
receive sides. 

In section II, we first recall the optimal power allocation for 
two-sided coordination vector channels (i.e. full vectoring) with self 
crosstalk and external noise under a total power constraint. The 
primal MIMO capacity optimization problem subject to a total power 
constraint coupled over the tones is transformed into a collection of 
per-tone unconstrained optimization problems using a dual problem 
formulation. We derive optimal transmitter and receiver structures 
(precoders and equalizers) in combination with power allocation 
which achieve MIMO channel capacity. Secondly, we devise the 
optimal power allocation for two-sided coordination vector channels 
with self crosstalk and external noise under per-modem total power 
constraints and spectral mask constraints, leading to a generalized 
SVD-based transmission (section III). Similar derivations are given 
for the special (or simpler) per-modem total power constraints case. 
Simulation results are given for bonded VDSL2 systems with external 
noise coming from ADSL2+ or VDSL2 disturbing lines, along with 
a comparison between algorithms with one-sided signal coordination 
either only at the transmit side or the receive side (section IV). 
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II. Total power constraint 

In this paragraph, we recall the classical SVD-based algorithm 
with optimal power allocation for two-sided coordination vector 
channels with self crosstalk and external noise under a total power 
constraint. We assume that transmitters use Discrete Multi-Tone 
(DMT) modulation with a cyclic prefix longer than the maximum 
delay spread of the channel. As shown in [8], there is more correlation 
in the spatial domain than the frequency domain. In this paper, we 
exploit the noise correlation in the spatial domain and we assume that 
the external noise is synchronized with the MIMO system. Therefore, 
there is no correlation in the frequency domain and the external noise 
is decoupled over the tones. The transmission over one tone can then 
be modelled as: 



y, = H,x, 



: = l. 



(1) 



where N c is the number of subcarriers, Xi is the vector of N 
transmitted signals on tone i, yt the received signal vector, H; the 
N x N MIMO channel matrix and the vector of noise containing 
Additive White Gaussian Noise (AWGN) and external noise (alien 
crosstalk, radio frequency interference,...). The primal problem of 
finding optimal PSD's for a MIMO binder under a total power 
constraint P tot is: 



max C(*i)i=i...Ar c 

(*i)i=l...W c 

subject to E Troce(*i) < P*' 

i=l 

*« y 0, i = 1 . . . N c 



(2) 



with 3?i the covariance matrix of transmitted symbols 3>j = E{k z x^] 
over tone i for the MIMO binder and with the objective function being 
the MIMO capacity summed over the N c tones [9]: 

C(* i )<=i... J v c = X><?2 [det (i + H.SiHfR- 1 )] (3) 

i=l 

Here, Ri is the covariance matrix of the noise Ri = E[niiif]. The 
idea of dual decomposition is to solve (2) via its Lagrangian [10]. 
The Lagrangian decouples into a set of N c smaller problem, thus 
reducing the complexity of equation (2). The dual objective function 



F(A) = max £(A, (*i) i= i.. 

(*i)i=l...N c 



(4) 



with 



£(A,(*i)i=i...* ) 



E 

»=i 



log 2 [det (I + H^HfR- 1 )] 



-\Trace(3>i) + AP 



(5) 



with A the Lagrange multiplier. The dual optimization problem is 
minimize f (A) 
subject to A > 



(6) 



Because the dual function is convex in A, standard convex op- 
timization results guarantee that the primal problem (2) and the 
dual problem (6) have the same solution [11]. Indeed, the objective 
and constraint functions are differentiable and the Slater's conditions 
are satisfied, therefore the duality gap is zero and the minimum 
of the dual function corresponds to the global optimum of the 
primal problem [10]. The search for the optimal A in (6) involves 



evaluations of the dual objective function (4), i.e. maximizations of 
the Lagrangian, which is decoupled over the tones for the given A. 
By exploiting the Cholesky decomposition R, = LjL^, where Li is 
a lower triangular matrix (whose inverse will be used to whiten the 
noise at the receive side), we obtain the following equation (using 
the property det (I + AB) = det (I + BA)): 



C(X, (*i) i= i...jv ) = E log 2 [detfl + L^H^HfL-r")] 

»=i 



-\Trace(4>i) + AP 



(7) 

The Singular Value Decomposition (SVD) of the withened channel 
L s _1 Hi = UiDiV^ transforms the inital channel into a product 
between two unitary matrices U», Vf and a diagonal matrix of 
singular values D;. 



C(X, (*<)i=i...jv c ) = E l°92 [det (I + UiDiVf * i V i D j Uf )] 

-\Trace(&i) + \P tot 
By setting <&i = Vf^jV, we can rewrite this as: 



£(A,(*i)i=i...JVc) = E [log* [det (l + D?*i)] 



-ATroce(*i) + AP 



(9) 



The off-diagonal elements in merely reduce the determinant 
owing to diagonal matrices 's and the property det(I + A) < 
Trace(l + A). Hence the optimal 4h is diagonal. In order to find 
the maximum, we compute the derivative of the function: 



d£(A,(*i) 



= ) 



d$, 



ln{2) 

The optimal is given by: 



diag 



(D 4 - a +* 4 )' 



I 



ln(2)X 



D" 



AI = 

(10) 

(ID 



where the [.]+ operation is inserted in order to obtain positive 
semi-definite *i's in formula (2). This is the well-known closed 
form waterfilling solution for MIMO systems 1 . The optimal power 
allocation consists of finding the optimal Lagrange multiplier which 
meets the total power constraint according to (11). The complete 
algorithm for the optimal power allocation under a total power 
constraint is given in the Annex A. The optimal Tx structure is given 
by the precoding matrix Vj while the optimal Rx structure is given by 



the equalizer matrix U f 
by: 



leading to parallel SISO systems as defined 



Uf L; Vi = DiXj + Uf L 4 1 n i 



(12) 



For practical implementations, we introduce the SNR gap T referred as the 
code gap in [12] which is the SNR multiplier required to achieve the target 
probability of error at the desired data rate. Considering the same T for the 
different virtual channels, the optimal 3>i is given by: 



Vi 



ln{2)X 



vf 
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with the optimal power allocation under a total power constraint 
driven by (11). 

III. PER-MODEM TOTAL POWER CONSTRAINTS AND SPECTRAL 

MASK CONSTRAINTS 
In the xDSL context, it is more relevant to consider a constraint 
on the power of each modem separately instead of a constraint on 
the power for all modems together. DSL standardization often defines 
spectral masks that each transmitter has to satisfy as well as the total 
power that each transmitter can transmit. In this section, we devise the 
optimal power allocation for two-sided coordination vector channels 
with self crosstalk and external noise under per-modem total power 
constraints (i.e. a single total power constraint for all tones per line) 
and spectral mask constraints as well as their corresponding optimal 
Tx/Rx structures. 

A. Optimal power allocation 

The primal problem of finding optimal PSD's for a MIMO binder 
under per-modem total power constraints Pj ot and spectral mask 
constraints is: 

max C(*i)i=i...jv c 

(*i)i = l...JV c 



subject to E < if* Vj 

*i y o, i = 1 . . . iV c 



(13) 



with the objective function being the MIMO capacity summed over 
the N c tones given by (3). Again, we can apply the idea of dual 
decomposition by decoupling the primal problem into N c smaller 
problems [10] considering the per-modem total power constraints and 
the spectral amsk constraints. The dual objective function is: 



F(A,Ai,...,Ajv c ) = max £(A,Ai , Ajv c , (*i)i=i...jv c ) 

(*t)t=l...iV c 

(14) 

with 



£(A, Ai, . . . , A Nc , (*<)»=i...jv c ) = 
E \log 2 [det (I + Hi*iHf Rr 1 )] 



A* 

I 

i=l 



(15) 



-Trace((A + A I )* I )J + Trace (Adiag(P* ot )) 

+ E Trace (k l diag(<j>™ aak - 3 ) 
i=i ^ 

The Lagrange multipliers corresponding to the per-modem total 
power constraints are contained in the diagonal matrix A = 
diag(Xi, . . . , Xn), the Lagrange multipliers corresponding to the 
spectral mask constraints for tone i are contained in the di- 
agonal matrix Ai = diag(Xa, . . . , Aijv)- The diagonal matrix 
diag(<j>™ ask ' j ) = diag{<f>™ ask ' 1 ,...,<l>™ asKN ) corresponds to the 
spectral mask for user j and tone i. The dual optimization problem 
is: 



minimize 



F(A,Ai,...,Aji 



subject to [A]jj, [Ai]jj, . . . , [A.N ]jj > V? 



(16) 



The dual function is convex in A, Ai, . . . , Ajv c , therefore standard 
convex optimization results guarantee that the primal problem (13) 
and the dual problem (16) have the same solution [11]. The La- 
grangian is differentiable and the Slater's conditions are satisfied, 
therefore the duality gap is zero and the minimum of the dual 



function corresponds to the global optimum of the primal problem 
[10]. The search for the optimal A, Ai, . . . , Ajv c in (16) involves 
evaluations of the dual objective function (14), i.e. maximizations 
of the Lagrangian, which is decoupled over the tones for the given 
matrices A, Ai , . . . , Ajv c . By exploiting the Cholesky decomposition 
Ri = LjLf , where Li is a lower triangular matrix (whose inverse 
will be used to whiten the noise at the receive side), we obtain the 
following equation (using the property det(l+ AB) = efet(I+BA)): 



£(A, Ai,. . . , Ajv c , (*i) i= i...jv c ) = 

E f logi [det (I + Lr 1 H i * i HfLr»)] 

i=i 



(17) 



-Trace{{\ + A,)*i)J + Trace (&diag(Pj ot )) 
+ E Trace (A l diag(^ ask - 3 )) 

i = l V ' 

The Singular Value Decomposition (SVD) of the whitened channel 
Lr 1 H i (A + Ai)" 1/2 = UiDiVf transforms the initial channel 
into a product between two unitary matrices Ui, ~Vf and a diagonal 
matrix of singular values Di. 

£(A, Ai,. . . , A Nc , (*i) i= i...jv c ) = 



E No S2 [det(I + UiDiVf (A + A0 1/2 *i 



(A + A i ) 1 / 2 V l D I Uf )] - Trace(Vf(A + Ai) 1 / 2 *, 



(18) 



(A + A0 1/2 V 8 )J + Trace (Adiag{Pj ot )) 
+ E \ Trace (kidiag{<l>™ ask ' j )) 



By setting # ; = Vf (A+Ai) 1/2 *i(A+A l ) 1/2 V I we can rewrite 
(18) as: 



£(A, Ai,..., Ajv c , (*i) i= i.. ,jv c ) = E yogi [det (^I + D 2 *^] 
-Trace{&i)} + Trace (Adiag(Pj ot )) 

(19) 

In order to find the maximum, we compute the derivative of the 
Lagrangian: 



- E Trace {A l diag{<j>™ ask ' i )) 
i=l ^ ' 



dC(A, Ai, . . . , Ajv c , (*i)i=i...jv c ) 



d® l 



D: 



(20) 



1 = 



The optimal is given by: 



*i = (A + Ai) 



-1/2 



Vi 



I 



ln(2) 



DJ 



Vf(A + A,)- 1/2 



(21) 



where the [.] + operation is inserted in order to obtain positive semi- 
definite 3>i's in formula (13). One can note that the precoder formulas 
are a function of the Lagrange multipliers A, Ai, . . . , Ajv c 's. This 
is the new generalized SVD-based closed form solution for MIMO 
systems under per-modem total power constraints and spectral mask 
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constraints 



B. Optimal Tx/Rx structure 

The Tx/Rx structure is obtained as follows. The first step is to 
find the optimal Lagrange multipliers defined for the per-modem 
total power constraint and the spectral mask constraints according 
to the dual objective function F(A, Ai, . . . , Ajv c ). As the function 
is continuous differentiable, the search algorithm can use a gradient- 
descent like method to find the optimal Lagrange multipliers and is 
guaranteed to converge. The algorithm tries to converge under the 
per-modem total power constraints over the tones and inside this 
optimization tries to converge on a per-tone basis to also satisfy 
the spectral mask constraints. The complete algorithm of power 
allocation under per-modem total power constraints and spectral mask 
constraints is given in the Annex B. After calculating the optimal 
Lagrange multipliers, we can calculate for each tone the SVD of the 
whitened channel scaled by the Lagrange multipliers L~ 1 H,(A opt -|- 
Ai,opt) -1 ' 2 = UiD;Vf and multiply the transmitted symbols by 
(A opt + Ai, opt ) _1 ' 2 Vi and the received symbols by Uf leading to: 

Uf f L~ 1 y i = Uf 1 H i (A opt + A i , op( )" 1/2 V i x i + Uf L" 1 ^ 

(22) 

This leads to parallel SISO systems as defined by: 

Uf L~ 1 y ! = D 8 x ; + Uf L-V (23) 

with the optimal power allocation under per-modem total power 
constraints and spectral mask constraints driven by (21). 

C. Per-modem total power constraints 

In this paragraph we recall the main steps for finding the optimal 
power allocation and optimal Tx/Rx structures with self crosstalk and 
external noise under per-modem total power constraints. We give the 
SVD-based algorithm under per-modem total power constraints for 
the readers who are interested in the derivations without referring to 
the more complicated SVD-based algorithm under per-modem total 
power constraints and spectral mask constraints 3 . The primal problem 
of finding optimal PSD's for a MIMO binder under per-modem total 
power constraints is: 

max C(*j)i=i...jv 

(*i)t = l...W 

subject to E < Pj ' W 

i=l 

t 0, i = 1 ■ ■ ■ N c 
The dual objective function is: 



F(A) = max £(A, (* i ) i =i... JVc ) 

(*i)>=l...N c 



(24) 



(25) 



with A a diagonal matrix of Lagrange multipliers 

diag(\ 1 , . . . , Xn) and 

2 For practical implementations, we again introduce the SNR gap T referred 
as the code gap in [12] which is the SNR multiplier required to achieve the 
target probability of error at the desired data rate. Considering the same T 
for the different virtual channels, the optimal &i is given by: 



*i = (A + A,) 1/2 V, 



I 

Zn(2) 



TD7 



Vf(A + Ai) 



-1/2 



3 The optimal power allocation under per-modem total power constraints 
can be found directly from (21) by setting Ai, . . . , Ajv c = 0. Moreover, the 
optimal power allocation under a total power constraint can be found directly 
from (21) by setting Ai, . . . , Ajv =0 and A = AI 



£(A, (*i)i=i...jv c ) = E log* [de^I + H^HfR- 1 )] 
-Trace(A$i) + Trace (Adiag(P<j ot )) 



The dual optimization problem is: 



minimize 

A 



F(A) 



subject to Aj > Vj 



(26) 



(27) 



The search for the optimal A involves evaluations of the dual 
objective function, i.e. maximizations of the Lagrangian, which 
is decoupled over the tones for a given set A/s. By exploiting 
the Cholesky decomposition Ri = L»Lf , by defining the (A- 
dependent) SVD L 8 _1 HiA" 1/2 = UjDjVf and by setting = 
Vf A 1 ^ 2 #iA 1/ ' 2 Vi, we can reformulate the optimization problem 
as: 



C(A, (h)i=i...N c ) = E (log, [det (i + D 2 *,)] 
-Trace(^i) ] + Trace {Adiag{P t J ot )) 



(28) 



We compute the derivative of the function in order to find the 
maximum, therefore the optimal power allocation is given by: 



-1/2 



I 



ln{2) 



DJ 



Vf A 



-1/2 



(29) 



The complete algorithm of power allocation under per-modem total 
power constraints is given in the Annex C. After calculating the 
optimal Lagrange multipliers, we can calculate for each tone the 
SVD of the whitened channel scaled by the Lagrange multipliers 
L~ 1 HiA~ p 1 / 2 = UjDjVf , where A opt is the optimal setting for 
the Lagrange multipliers, we multiply the transmitted symbols by 
A~p/ 2 Vi and the received symbols by Uf leading to: 



Uf L'Vi = Uf L~ 1 H i A~ p 1 / 2 ViX i + Uf L-V (30) 
This leads to parallel SISO systems as defined by: 

UfLrVi = D iXl + Uf L^nj (31) 

with the optimal power allocation under per-modem total power 
constraints driven by (29). 

D. Remark 

The power allocation problem under a total power constraint and 
spectral mask constraints can be solved similarly to the problem 
of power allocation under per-modem total power constraints and 
spectral mask constraints. In this case A = AI. The same derivation 
as in the previous section can be given. The optimal power allocation 
under total power constraint is then given by: 



(AI-r-Ai)- 1/2 V 8 



I 



ln{2) 



D" 



Vf (AI-r-Ai) _1/2 (32) 



The PSD's formula and the algorithm description can be easily 
modified accordingly. 
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Fig. 1. Rates obtained in a downlink scenario of the SVD-based algorithm 
with per-modem total power constraint and the Diagonalizing Precoder for 
ADSL2+ and VDSL2 disturbing lines 
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Fig. 2. Rates obtained in an uplink scenario of the SVD-based algorithm 
with per-modem total power constraint and the ZF receiver for ADSL2+ and 
VDSL2 disturbing lines 



IV. Results 

Far End Crosstalk (FEXT) and Near End Crosstalk (NEXT) models 
are well known in the literature for xDSL [7]. The NEXT and FEXT 
models for ADSL/HDSL may be no longer applicable to VDSL2 
due to much larger bandwidth. In this paper, the simulations results 
are obtained on measured channels from a France Telecom binder 
with 8 lines of 800 meters and external noise coming from 400 
meters lines. We look at the performance of VDSL2 lines with two 
different sources of external noise. The two sources of external noise 
consist of ADSL2+ or VDSL2 disturbing lines having their PSD's at 
spectral masks [13], [14]. This gives an indication of the statistical 
variations of the achievable bit-rates of the presented algorithms 
because different combinations of pairs may be bundled. Recently, 
crosstalk channels for VDSL2 have been characterized as parametric 
models [15]. Therefore, the optimal Tx/Rx structures and power 
allocation algorithms proposed in this article can also be used in 
such generic models 4 . 

We use spectral masks for VDSL2 Fiber To The exchange (FTTex) 
as described in [14], with SNR gap r=10.8 dB (Shannon gap=9.8 
dB, margin=6 dB and coding gain=5 dB) to achieve the target BER, 
an AWGN of -140 dBm/Hz and maximum transmit power P?°*=14.5 
dBm per line. The power spectrum of the disturbing system (ADSL2+ 
or VDSL2) is set to its spectral mask [13], [14]. The frequency 
range is from to 12 MHz with 4.3125 kHz spacing between 
subcarriers and 4 kHz symbol rate. The FDD band plan of VDSL2 
corresponds to 2 frequency bands in the downlink scenario which are 
138kHz-3.75MHz and 5.2MHz-8.5MHz. In the uplink scenario, this 
corresponds to 3 frequency bands 25kHz-138kHz, 3.75MHz-5.2MHz 
and 8.5MHz-12MHz. The processing of the bundled systems is the 
same for both cases. 

Fig.l shows the comparison between the two-sided coordination 
vector channels SVD-based algorithm and the one-sided coordination 
BC Diagonalizing Precoder (DP) algorithm with per-modem total 
power constraint in a downlink scenario [4]. The BC DP corresponds 
to a scaled version of the ZF precoder by the diagonal elements of 

4 Note that from the implementation point of view, once the optimal covari- 
ance matrices (^i)i— i...jv c are determined, the transmitted data symbols 
can be constructed as follows: 

1) The N X 1 vector of the M-QAM data symbols is precoded using 
the N X N Ci matrix, i.e. Xj = C;s;, such that £[x;xf] = *, = C,Cf 
from Cholesky decomposition (as 3>i is a positive semi-definite matrix). 

2) Then x x will be sent on the N lines (the j th element of x^ will be sent 
on the j* line). 



the channel matrix. The length of the bonded lines are 800 meters 
while the coupling between the bonded lines and the disturbing 
lines occurs the last 400 meters to the Customer Premise Equipment 
(CPE). ADSL2+ disturbing lines and VDSL2 disturbing lines whose 
PSD's are set to their respective spectral masks are simulated. A 
MIMO binder of 8 lines is used, with the number of coordinated 
pairs going from 1 to 8 and the number of disturbing lines from 
7 to respectively. In each case, all (^) combinations are used 
to provide an average bit-rate. When the number of coordinated 
pairs equals to 8, there is no external noise and thus this pro- 
vides maximum performance. The SVD-based transmission with per- 
modem total power constraint performs better than the DP with per- 
modem constraint owing to the exploitation of the noise correlation 
by whitening. Moreover, the higher the number of coordinated pairs, 
the higher the improvement in terms of bit-rate. In fact, the whitening 
process provides more cancellation performance of the disturbing 
lines when the number of coordinated pairs is higher than the number 
of disturbing lines. Although the spectral mask of ADSL2+ is much 
higher than VDSL2 for frequencies up to 2.2 MHz, the VDSL2 
disturbing lines have much more impact on the considered bonded 
lines because the crosstalk increases as the frequency increases. 

We have simulated the one-sided coordination BC Minimum 
Mean Square Error-Dirty Paper Coding (MMSE-DPC) algorithm in a 
downlink scenario with optimal power allocation found by exhaustive 
search without external noise cancellation since noise whitening is not 
possible at the transmit side [2]. For the 2 user case the (*) channels 
have been processed and they lead to an average bit-rate of 14.92 
Mbps in the case of VDSL2 disturbing lines compared to 15.60 Mbps 
for the SVD-based solution with per-modem total power constraints. 
Therefore the external noise decreases the performance of the one- 
sided BC MMSE-DPC algorithm in a downlink scenario. The SVD- 
based algorithm provides an upper bound on the achievable capacity, 
thus it exploits the equivalent whitened channel matrix where optimal 
power allocation is found by a closed form formula contrary to one- 
sided coordination BC MMSE-DPC where optimal power allocation 
is found by exhaustive search. 

Fig.2 shows the comparison between the two-sided coordination 
vector channels SVD-based algorithm and the one-sided coordination 
MAC ZF receiver with per-modem total power constraint in an uplink 
scenario [3]. Similar conclusions can be told for this comparison, 
however even if the whitening step is possible at the receive side, the 
ZF equalizer can't take advantage of the equivalent channel as the 
SVD scheme does. Indeed, the ZF structure consists of an inversion of 



6 



IEEE TRANSACTIONS ON COMMUNICATIONS, VOL. XX, NO. Y, MONTH 2007 



— ADSL2+ average 
-HADSL2+ minimum 
-*-ADSL2+ maximum 
---VDSL2 average 
-+-VDSL2 minimum 
-*- VDSL2 maximum 






4 ■ — 



1 2345678 
Number coordinated lines 

Fig. 3. Rates obtained in a downlink scenario for the different permutations 
of SVD algorithm under per-modem total power constraints and spectral mask 
constraints for ADSL2+ and VDSL2 disturbing lines 
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Fig. 6. Rates obtained in an uplink scenario of the different SVD-based 
algorithms under total power constraint, per-modem total power constraints 
and spectral mask constraints and per-modem total power constraints using (a 
posteriori) truncation by the spectral mask for ADSL2+ and VDSL2 disturbing 
lines 
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Fig. 4. Rates obtained in an uplink scenario for the different permutations of 
SVD algorithm under per-modem total power constraints and spectral mask 
constraints for ADSL2+ and VDSL2 disturbing lines 
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Fig. 5. Rates obtained in a downlink scenario of the different SVD-based 
algorithms under total power constraint, per-modem total power constraints 
and spectral mask constraints and per-modem total power constraints using (a 
posteriori) truncation by the spectral mask for ADSL2+ and VDSL2 disturbing 
lines 



the channel matrix, and the performance is not changed by inverting 
the whitened channel. Moreover, there is a very small difference 
between ZF and SVD schemes for ADSL2+ crosstalkers because 
they don't harm the VDSL2 lines. Indeed, the bandwidth involved 
in the uplink scenario is the 25kHz- 13 8kHz bandwidth, therefore 
SISO schemes with 7 ADSL2+ crosstalkers and 8x8 MIMO schemes 
with no ADSL2+ crosstalkers have almost the same performance. 
Contrary to ADSL2+ crosstalkers, VDSL2 crosstalkers decrease the 
performance of the VDSL2 lines. For the 2 user case, the one-sided 
coordination MAC MMSE-Decision Feedback Equalizer (MMSE- 
DFE) algorithm leads to an average bit-rate of 4.20 Mbps in the case 
of VDSL2 disturbing lines, which is similar to the average bit-rate of 
the SVD-based algorithm. In this case, the whitening is possible at the 
receive side and the MMSE-DFE receiver can exploit the equivalent 
channel. 

Fig.3 and Fig.4 show the bit-rate performance of the two-sided 
coordination vector channels SVD algorithm under per-modem total 
power constraints and spectral mask constraints for the downlink and 
the uplink scenario with ADSL2+ and VDSL2 disturbing lines. The 
(^) combinations are used to provide an average bit-rate, a minimum 
bit-rate and a maximum bit-rate, thus the variations in bit-rates using 
all different permutations. One can observe that the range of bit- 
rates diminishes as the number of coordinated lines increases. The 
comments given in Fig.l and Fig.2 also apply to these figures. 

Fig.5 and Fig.6 give the bit-rate performance of the two-sided co- 
ordination vector channels SVD-based algorithms under a total power 
constraint, per-modem total power constraints using (a posteriori) 
truncation by the spectral mask 5 , per-modem total power constraints 
and spectral mask constraints for the downlink and the uplink scenario 
with ADSL2+ and VDSL2 disturbing lines. A total power constraint 
gives an extra degree of freedom compared to per-modem total power 
constraints and thus achieves better performance. Similar comments 
from the previous figures could be told. We observe that there 
are very small differences between bit-rates under a total power 
constraint and bit-rates under per-modem total power constraints. For 
the SVD algorithm with per-modem total power constraints using (a 
posteriori) truncation by the spectral mask, the spectral masks are 

5 The truncation by the spectral mask cuts the optimal PSD's found by the 
SVD-based algorithm under per-modem total power constraints and do not 
distribute to further tones the power loss due to the truncation when the PSD's 
are higher than the spectral mask 
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directly applied to the optimal PSD's under per-modem total power 
constraints. There is a noticeable difference between the performance 
of the SVD algorithm using truncation by spectral mask and the SVD 
algorithm under spectral mask constraints especially in the uplink 
scenario. In fact, this optimization process allows a better distribution 
of the power over the tones by setting the optimal PSD's larger than 
the spectral mask at the spectral mask and thus saving power for 
other tones. 

We can notice that the gain of the proposed techniques originates 
from the MIMO diversity gain and does not originates from water- 
filling on the considered loops. In fact, the eigenvalues used in the 
SVD-based algorithms do not have a significant impact on the optimal 
power allocation. The inverse of eigenvalues is in the order of 1CT 8 
and the water level is around 10" 6 for the considered loops and the 
total power constraint in (11) (the same behaviour can be observed 
under per-modem total power constraints and mask constraints in 
(21)). As the heart of the waterfilling solution reduces the water level 
by the inverse of the eigenvalues, we observe only small variations 
around the water level and the waterfilling operation leads to flat 
PSD's. Destroying this diagonal dominance by increasing artificially 
the crosstalk before the withening operation, results in a capacity gain 
originating both from MIMO diversity and the waterfilling gains. The 
resulted optimal power allocation (not included in the manuscript for 
the space limitations) is no longer flat for the two lines. This suggests 
that these SVD-based algorithms could have a significant impact 
in the wireless context where the channel matrix is not diagonally 
dominant. 

In general, in VDSL2 scenario's with a binder with equal length 
cables, we can expect that per-modem total power constraints and 
spectral mask constraints do not degrade performance, while these 
constraints do lead to more practical, implementable and standards- 
compliant solutions. Moreover, power allocation algorithms under 
spectral mask constraints can provide better performance compared 
to a simpler power allocation procedure where the spectral mask 
constraints are first removed from the optimization problem, and then 
imposed onto the computed PSD's. This is especially so in uplink 
scenario's or even more in a scenario where optimal PSD's could be 
much larger than spectral mask PSD's. 

V. CONCLUSION 
In this paper, we have investigated the problem of finding the 
optimal power allocation in MIMO xDSL systems under self crosstalk 
and external noise and with two-sided coordination (full vectoring). 
Optimal Tx/Rx structures and power allocation algorithms have been 
devised under practical limitations from xDSL systems, namely per- 
modem total power constraints and spectral mask constraints, leading 
to a generalized SVD-based transmission. Simulation results were 
given for bonded VDSL2 systems with external noise coming from 
ADSL2+ or VDSL2 disturbing lines, along with a comparison with 
algorithms with one-sided signal coordination, either only at the trans- 
mit side or the receive side. The two-sided coordination SVD-based 
algorithms then provide a performance upper bound for the existing 
one-sided coordination MAC ZF, BC DP, MAC MMSE-DFE or BC 
MMSE-DPC algorithms under self crosstalk and external noise. The 
simulation results also showed that adding per-modem total power 
constraints and spectral mask constraints did not significantly reduce 
the bit-rate compared to the case where only a total power constraint 
is imposed owing to the SVD-based transmission. The optimal power 
allocation algorithms under spectral mask constraints also provides 
better performance compared to a simpler power allocation procedure 
where the spectral mask constraints are first removed from the 
optimization problem, and then imposed onto the computed PSD's. 
An extension of this work could be aimed at canceling external 



noise using the noise correlation between differential-mode and 
common-mode signals, or between different tones when considering 
asynchronous external noise. 



Annex 



This Annex provide the algorithms for finding the optimal La- 
grange multipliers for the two-sided coordination vector channels 
SVD-based algorithms under a total power constraint, per-modem 
total power constraints and spectral mask constraints and per-modem 
total power constraints respectively. 



A. Total power constraint 

The following algorithm provides the optimal power allocation for 
the two-sided coordination vector channels under self crosstalk and 
external noise using an SVD-based algorithm. The later tries to find 
in an iterative way the optimal Lagrange multiplier to meet the total 
power constraint. 



Algorithm 1 Total power constraint 
init A = 1 
init step = 2 
init 6 = 

while | Trace(&i) - P tot \ > tolerance 

i=l 

if J2 Trace(&i) - P tot < 

i=l 

6 = 6+1 

A = X/step 

step = step - l/2 b 

end if 

A = A * step 

*»= V »[OT- D ^] + V ^ Vi 
end while 



B. Per-modem total power constraints and spectral mask constraints 

The following algorithm provides the optimal power allocation for 
the two-sided coordination vector channels under self crosstalk and 
external noise using an SVD-based algorithm. The later tries to find in 
an iterative way the optimal Lagrange multipliers to meet per-modem 
total power constraints and spectral mask constraints. 
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+ Vf (A + Ai) 1/2 Vi 



Algorithm 2 Per-modem total power constraints and spectral mask 

constraints 

init Xj = 1 Vj 
init stepj = 2 Vj 
init bj = Vj 

init *i = (A + Ai) -1/,2 Vi [j^-D- 2 

while | ^ - -P, tot | < tolerance 3j 

'"for i=l to W c 

init Xj = 1 Vj 
init step,,- = 2 Vj 
init 6j = V j 
while < iterations 
for j=l to N 



(A + A.r^V, - D- 2 ] + Vf (A + A,)" 

if [*;]„■ > 



1/2 



Vi 



Vi 



[*«]. 



Vi 



end if 

if - ^° 3fc < Vi 



= Xj/stepj 
stepj = step 3 - 

end if 

Xj = Aj * stepj 

end for 
end while 

end for 
for j=l to N 

if E " P) ot < 

i=i 

6^ =bj + l 

Xj — Xj/stepj 

stepj = stepj — l/2 b i 

end if 

Xj = Aj * sfepj 

end for 



■ 1/2 6 



= (A + A,)' 



Vi 

end while 



1 



(2) 



Vf(A + Ai) 



-1/2 



C. Per-modem total power constraints 



The following algorithm provides the optimal power allocation for 
the two-sided coordination vector channels under self crosstalk and 
external noise using an SVD-based algorithm. The later tries to find in 
an iterative way the optimal Lagrange multipliers to meet per-modem 
total power constraints. 



/2 



1U lOl 

>=A- 1/2 V 8 [^-D- 2 ] VfA- 



/2 



Algorithm 3 Per-modem total power constraints 
init X-j = 1 Vj 
init stepj = 2 Vj 
init bj = Vj 

init ^A-^v^-D-J + VfA" 1 

while | J2 [$i]jj - Pj ot \ > tolerance Bj 
'for j=l to N 

if 2[**]«--pr<° 

i=l 

6^ = foj + 1 

Aj = Xj/stepj 

stepj = stepj - l/2 b i 

end if 

Aj = Aj * , 

end for 
end while 
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